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We report on the spatio-temporal dynamics of the interface in shear-banding flow of a wormlike 
micellar system (cetyltrimethylammonium bromide and sodium nitrate in water) during a start-up 
experiment. Using the scattering properties of the induced structures, we demonstrate the existence 
of an instability of the interface between bands along the vorticity direction. Different regimes of 
spatio-temporal dynamics of the interface are indentified along the stress plateau. We build a model 
based on the flow symetry which qualitatively describes the observed patterns. 
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Many complex fluids of various microstructures 
strongly exhibit non linear properties under simple shear 
flows. When the characteristic times of the flow (~ I/7) 
and of the system are similar, the system is liable to 
undergo instabilities and flow- induced phase transitions. 
Such coupling between flow and structure has been ob- 
served for example, in wormlike micellar solutions 0, 
lyotropic lamellar phases |2[ , telechelicpolymer networks 
Pj, soft glassy materials Q or foams pj. 
For a system showing a shear-thinning behaviour, the 
mechanical signature, in the measured flow curve (a — 
f(j)), of this type of phenomena is the stress plateau 
separating low and high viscosity branches. This stress 
plateau is at the origin of a shear-banding transition from 
a flow at first homogeneous towards a stationary state 
where two layers supporting different shear rates coexist. 
A change of the applied shear rate only affects the rel- 
ative proportion of each layer, the width of the induced 
band increasing linearly with the macroscopic shear rate 



Among the complex fluids with stress plateau, the sys- 
tems of reversible giant micelles have been the object of 
the most intensive surveys (see for a review). Very 
recently, the use of dynamic light scattering to measure 
velocity profiles brought confirmation of the very sim- 
ple shear-banding scenario previously mentioned for the 
widely studied CPCl/NaSal solution Q. However, the 
authors mention the existence of temporal fluctuations 
in the highly sheared band. More complex pictures have 
also emerged with the development of time- and spatially- 
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resolved techniques such as rapid nuclear magnetic res- 
onance (NMR) velocimetry 8] and ultrasonic velocime- 
try 0, with for example, the observation of periodic or 
quasirandom oscillations of the interface position driven 
by wall slip. Moreover, the flow birefringence, technique 
sensitive to molecular alignment, revealed complex kinet- 
ics of banding formation and strongly differing banding 
organizations from one sample to another [lfj. 
Some of these behaviors are reproduced, at least quali- 
tatively, by several theoretical works, which predict fluc- 
tuating or chaotic flows, taking into account a coupling 
between flow and microstruture via the concentration or 
the micellar length [llj. They also motivated a recent 
theoretical work, in which the author, using the nonlocal 
Johnson-Segalman model, studies the stability of planar 
shear banded flow with respect to small perturbations 
with wavevectors in the interfacial plane 12]. An inter- 
facial instability is predicted in the plateau domain for 
almost all shear rates, the interface profile being then 
modulated by a wave with wavevectors confined along 
the flow direction. In the present letter, we report on the 
formation kinetics of the banding structure and follow 
particularly the spatio-temporal dynamics of the inter- 
face between bands for a sample of cetyltrimethylammo- 
nium bromide (CTAB) and sodium nitrate (NaN03)^j|. 
Direct visualizations of the gap of the Couette cell in the 
plane vorticity/ velocity gradient and variations of the 
shear stress are recorded simultaneouly during a start- 
up experiment. Using the scattering properties of the 
induced structures, we show that the interface between 
the two bands becomes unstable with wavevector in the 
vorticity direction. Different patterns of spatio-temporal 
dynamics of the interface are observed, depending on the 
applied shear rate. Using arguments based on the sym- 
metry of the flow, we show that the interface position is 
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FIG. 1: (a) Steady state flow curves obtained under con- 
trolled strain rate (o, 250s per data point) and controlled 
shear stress (A, 48s per data point). The dotted lines point 
out the bounded characteristic values of j and a associated 
with the stress plateau, (b) Variation of the shear stress 
versus time after the onset of flow for applied shear rates 
7 = 10s" 1 and j — 30 s . The inset shows the response on 
a shorter time scale. 

described by a generic amplitude equation of Kuramoto- 
Sivashinsky type, which reproduces qualitatively the ob- 
served patterns [lj, [l5| . 

The wormlike micellar system chosen in this work is 
made of CTAB at 11% wt. and NaN0 3 (0.405 M) in dis- 
tilled water at T=30°C. At this concentration, far from 
the transition concentration towards an orientated ne- 
matic phase at rest, the micelles form a highly entan- 
gled network. The linear viscoelastic behavior of the 
solution is Maxwellian with a plateau modulus Go = 
238 ± 7Pa and a terminal relaxation time tr — 0.15 ± 
0.02s. The rheological measurements are performed on 
a stress-controlled rheometer (Physica MCR500) work- 
ing in strain-controlled mode. A home-build transpar- 
ent Couette cell with a height of 40 mm is used for the 
experiments. The inner rotating cylinder has a radius 
i?i = 13.33 mm providing a gap thickness e = 1.13 mm. 
The Couette cell is designed in a way which allows the di- 
rect observation of the gap in the plane vorticity/velocity 
gradient. We form a thin laser sheet in that plane with a 
cylindrical lens and record the scattered intensity at 90° 
with a CCD camera. 

The steady-state rheological behavior of our sample in 
strain- and stress-controlled modes of the rheometer is 
drawn in the semi- logarithmic plot in Fig l.a. The flow 
curve is composed of two increasing branches separated 
by a stress plateau extending from ji = 7 ± 0.5s -1 to 
A fh = 140 ± 10s" 1 , characteristic of the emergence of 
heterogeneous flow of shear-banding type. This plateau 
presents a significant positive slope partly due to the non 
homogeneity of the Couettegeometry. In fact, the curva- 
ture effects, analysed as in y\ do not allow for a complete 
explaination of the slope. Concentration effects between 
shear bands could also make the plateau steeper [l6j . 
We shall now focus on the kinetics of formation of the 
heterogeneous flow. At t=0, a shear rate is applied to 
the sample initially at rest and the evolution of the shear 
stress and direct observations of the gap of the Couette 



cell in the plane vorticity/velocity gradient are recorded 
simultaneously as function of time until steady-state is 
achieved. Fig. l.b illustrates the variations with time 
of the shear stress for two shear rates associated with 
the plateau region (7 > 7;). The stress a(t) follows the 
generic behavior observed on various semi-dilute micel- 
lar systems |l0t Il7j undergoing shear-banding transition 
induced by the flow: an overshoot at short times, then a 
slow sigmoidal relaxation (7 = lOs^ 1 ) or damped oscil- 
lations (7 = 30s -1 )) and finally a small undershoot that 
continues the transient response towards the stationnary 
state. 

Let us correlate such an evolution with the direct visual- 
izations of the gap (Fig. 2). The sample at rest does not 
present any particular scattering properties as observed 
on the first photo (t — 0~) where the gap appears dark. 
Just after the inception of flow, the solution becomes 
turbid over the entire gap and strongly scatters the in- 
cident light in all directions. This strong increase of the 
turbidity is associated with the elastic overshoot in the 
shear stress response (see photo 1). The entangled mi- 
cellar network is strongly stretched leading to concentra- 
tion fluctuations. Then, the liquid becomes transparent 
again near the fixed wall of the Couette cell (photo 2), 
this phenomenom correponding to the rapid relaxation 
of the stress overshoot. Let us note however, that the 
turbidity does not relax all over the gap since a scatter- 
ing band persists against the rotating inner cylinder, the 
interface between bands being diffused (photo 4). The 
crucial point here, is that the induced structures are nu- 
cleated during the stress overshoot and in all the gap of 
the cell so as to the induced band is already formed just 
after the stress overshoot. The enhancement of the con- 
centration fluctuations acts as the precursor of the tran- 
sition. Both bands have differing optical properties: the 
induced band is very turbid, indicating that it contains 
structures of mesoscopic size which scatters the light in 
all directions. Such a difference allows us to perflectly dis- 
tinguish the interface between bands. Let us note that 
the turbidity is associated with birefringence banding. 
From t = 10s the interface becomes clearly sharp (photo 
5). This lasts for about t = 20s (photo 6) and then 
we observe a completely unexpected behavior: the sharp 
interface between the two bands begins to destabilize it- 
self along the vorticity axis. The instability grows up to 
t = 29s (photo 9), namely when the stress undershoot 
occurs in the transient response. Beyond, the interface 
seems to adopt a stationnary periodic profile with a well- 
defined characteristic wavelength. An extension of the 
field of observation shows that this ripple of the interface 
spreads over all the height of the inner cylinder. At the 
edges of the cylinder, the ripple presents a minimum of 
amplitude. 

This destabilization process of the interface occurs in 
all the plateau region of the flow curve. In a general 
way, the wavelength and the amplitude of the instabil- 
ity clearly increase with the applied shear rate, except 
when 7 approaches jh where they diminish again. Let 
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FIG. 2: Views of the gap in the plane velocity gradient/vorticity taken out from the recording of the scattered intensity at 
different instants during the formation of the induced band. The sample is subjected to a shear rate of 30s~ . The letters A 
and B indicate respectively the inner and outer cylinders of the Couette device. Due to the compromise between the spatial 
resolution and the size of the field of observation, this latter is limited to 1.5cm in height and centered at halfway of the cell. 



us note that above 77,, the gap is completely filled by 
the induced band. Using a simple procedure of image 
analysis, we detect the interface profile on each frame. 
We represent in gray levels, the amplitude of the inter- 
face as function of time and space coordinate along the 
vorticity axis as illustrated in Fig. 3. a and b., where we 
summarized the complete behavior of the interface on 
spatio-temporal diagrams for 7 = 10s _1 and 7 = 30s" 1 
respectively. In first approximation, we adjust the in- 
terface profile by a sinusoidal function in order to esti- 
mate the wavelength and the amplitude of the instabil- 
ity. The observed patterns are very different for both the 
applied shear rates. At 10s _1 , the deformation of the 
interface occurs around t = 20 s and the waves formed 
seem to propagate along the z-coordinate, alternatively 
towards the top and the bottom of the cylinder, inter- 
acting in a complex scenario. We find for the amplitude 
and the mean wavelength, respectively 0.007±0. 003mm 
and 0.4±0.1mm. The evolution at 7 — 30s -1 appears 
simplier since there is no propagation in this case. The 
interface keeps a spatially stable profile in course of time 
with an amplitude of 0.05±0.01mm and a wavelength of 
2.4±0.02mm namely, more than two times the gap width. 
We checked by increasing the duration of the experiment 
up to 40 minutes that there is no coarsening between 
the "domains" (black and white) at very long times. A 
closer inspection of the diagram reveals the existence of 
time-periodic oscillations. We can distinguish two differ- 
ent frequencies: the first one, in the dark zones, is due 
to the mechanical imperfection intrinsic to the rheome- 
ter ; the second one, which is visible in the clear zones, 
is lower and results from the superposition of the default 
of coaxiality of the cell with a " beat" of the minimum of 
amplitude of waves. In other words, the amplitude of the 
interface instability oscillates with time. This could be 
due, among others, to tridimensional flow or to destabi- 
lization of the interface in the velocity direction [l^ . 



Moreover, from the asymptotic behavior we are able to 
deduce the proportion of the induced "phase" cf>h- the 
integration of the interface profile gives 4>t = 0.29 ± 0.02. 
The same calculus for each photo from t = 10s gives 
the same relative proportion. According to the simplest 
shear-banding scenario described by Cates's theory [(| 
and recently observed by Salmon et al 7], the width 
of the induced band increases linearly with the applied 
shear rate. In our case, the proportion of the induced 
band computed from the lever rule seems to be underes- 
timated : for 7 = 10s _1 and 7 = 30s _1 , we find respec- 
tively <f> h = 0.02 ± 0.01 and <j> h = 0.18 ± 0.02, values to 
be compared to 0.04 ± 0.01 and 0.30 ± 0.02. Processing 
of ultrasonic data on this sample is in progress in order 
to properly correlate the optical response and the local 
structure of the flow. 
We would like to propose a simple model that cap- 
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FIG. 3: Spatio-temporal evolution of the interface amplitude 
for 7 = 10s" 1 (a) and 7 = 30s" 1 (b) on a time scale which 
provides both kinetics of destabilization of the interface and 
the asymptotic behavior. The z-coordinate gives the spatial 
position along the vorticity axis. A linear gray scale (in /xm) 
is used, the black and white regions being associated respec- 
tively with minima and maxima of the interface amplitude. 
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FIG. 4: Spatio-temporal evolution of the interface position 
p obtained by numerical simulation of the model with (a) 
H = 0.07 (b) n = 0.2. 

tures the asymptotic dynamics of the fluid flow. As 
stated previously, when the fluid is forced with a shear 
rate that is inside the plateau region, the system spon- 
taneously coarse-graines into two zones separating two 
different shear rates. Here, as in lubrication approxima- 
tion for newtonian fluids, gradients with respect of the 
transversal coordinate y are higher than those with z the 
coordinate tangential to the cylinder axis. This length 
scales difference yields that the dynamics with respect 
to z coordinate are slow, and will consequently drive the 
interface position for asymptotic times. Let us note p the 
position of the interface separating these two zones. Tem- 
poral evolution of the front must reflect the symmetry of 
the experiment z — > —z, so only terms including spatial 
derivatives like d zz p, d zzzz p and {d z p) 2 are allowed at 
lowest order. The system is not invariant with respect 
to translation in the y direction, and the equation must 
present the broken symmetry : p —> p + constant. The 
asymptotic dynamics of the interface are given by: 

d t p = - p) - d zz p - d zzzz p + (d z p f . (1) 

Equation^ is written in non dimensional variables, and 
introduces the control parameter fj, in which all the phys- 
ical attributes of the experiment have been set. The first 
term of the r.h.s of the equation ^ just states that if the 
interface has a stationary location. Since the fluid is set 
into negative viscosity inside the plateau, the second term 
is a diffusion term with a negative sign. The next term, 
d zzzz p is a linear saturating term. The last term is the 



first non linear term that respects the symmetry at that 
order. This equation is known as the damped Kuramoto- 
Sivashinsky equation ^^|. When fi > 1/4 the solution 
p = 1 is a linearly stable solution. For \i < 1/4, the per- 
turbed interface becomes oscillatory in space, with a well 
defined wave-length. Finally, for fj, < 0.078, the dynamics 
are chaotic [lj| , see for example Fig. |U(a). When fi = 0, 
equation ^ reduces to the usual Kuramoto-Sivashinky 
equation, and exhibits complex spatio-temporal behav- 
ior that is used to describe phase turbulence 01 . The 
existence of stationary periodic pattern observed in the 
experiment is related to the very low aspect ratio of the 
Couette cell. Let us finally note that, in the model, when 
the interface becomes unstable, its position is systemat- 
ically shifted due to the non- linear term (d x p) 2 . This 
term has the same effect as the curvature of the Couette 
geometry. This means that the instability of the interface 
could also explain the "additional" slope of the plateau. 
In summary, our rheo-optical study of the interfacial dy- 
namics in shear-banding flow reveals the existence of a 
destabilization process of the interface between bands 
with wavevector in the vorticity direction. This observa- 
tion suggests that the flow field is three dimensional with 
the presence of recirculations. We show that there are 
different spatio-temporal dynamics of the interface along 
the stress plateau, with chaotic events at low shear rates, 
and spatially stable oscillatory behavior at higher strain 
rates. The physical origin of the destabilization process 
is still to be determined ; however, it could partly ex- 
plain some fluctuating behaviors highlighted recently in 
shear banded flows 0, H, an d it emphasizes the ne- 
cessity for a complete spatio-temporal description of the 
flow and the microstructural organization in the differ- 
ent planes. Further, we buid a model equation ^ based 
upon symctry arguments which described qualitatively 
the observed complex patterns. This model equation can 
also be derived using perturbation theory. To this end, 
it is necessary to adopt a model for the rheology of the 
fluid Ha. l20l . It is a computation that we did using the 
model|20| and we will report it elsewhere [2l| with the 
complete exploration of the stress plateau. 
The authors are grateful to G. Gregoire, O. Grefner, O. 
Cardoso, J.F. Berret, and J.L. Counord for the develop- 
ment of the transparent Couette cell. 
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